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Abstract 

Phase behavior of AB-multiblock copolymer melts which consists of chains with 
Laser distribution of A and B blocks have been investigated in the framework of the 
mean-field theory, where the polydispersity of copolymer is a function of two param- 
eters K and M. The influence of the Laser distribution on higher order correlation 
functions ( up to sixth order ) are computed for various values of K and M, and their 
contributions on the phase diagrams and phase coexistence are presented. It is shown 
that, with increasing polydispersity ( decreasing K and increasing M ) the transition 
lines of all phases shift upwards, consequently polydispersity destabilize the system. 
Keywordsrphase coexistence, Laser distribution, polydispersity, block length, 
Phase diagram, multiblock copolymer melts, correlation functions . 
PACs numbers:05.45.Ra, 05.45.Jn, 05.45.Tp 
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Recently several publications have been appeared in connection with block copolymer sys- 
tems from both theoretical ^, || and experimental @, |5|, § aspects. In AB multiblock 
copolymer melts, the immiscibility between A and B blocks induces self-assembly into vari- 
ous ordered microstructures. With decreasing the temperature, multiblock copolymers with 
repulsive interaction between different monomers undergo a microphase separation transition 
manifesting an energetic preference for contacts between similar monomers. In this work , 
the system will be studied in the weak segregation regime, where the separation between 
the A and B monomers is not complete. The theoretical investigations demonstrate that it 
is very possible that these (microphases separated) single phase are metastable and the free 
energy can be lowered even more by splitting the system into coexistence phase f?]. The 
properties of copolymer depend not only on the nature of the comonomers and the overall 
compositions, but also on the distribution of monomer units along the chain. As an example 
one can consider N-vinyl pyrrolidone (VP) and vinyl acetate (VA) copolymer. This copoly- 
mer can be synthesized in different ways resulting in the different homogeneity of backbone 
monomer distribution. Experimental studies show that the monomer unit distribution has 
effect on solution properties, particularly when the mole ratio VP to VA in the copolymer 
is close to unity, sequence distribution plays an important role in solution behavior ||. In 
this work as in Reference ||, we use the most comprehensive distribution, namely Laser 
distribution with two parameters K and M as the controller of degree of polydispersity. 
Polydispersity is an important parameter which can effect many physical properties of a 
polymer. Recently in an experimental work, the polydispersity lower than 1.0 is reported 
[|T0f . Also the theoretical models with polydispersity values of lower than unity can be useful 



for molecular weight distributions based on complex numbers ["J I 



The aim of the present paper is to investigate the effect of the parameters of the Laser distri- 
bution on the position of the phase transition lines, the stability of the various microstructre 
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and the position of the phase coexistence in multiblock copolymer. 

We consider an incompressible multiblock copolymer melt in which both the lengths of 
A-blocks and B-blocks satisfy Laser distribution. The phase behavior of the multiblock 
copolymer melt is calculated in the mean-field theory. By using this theory it is shown that 
there are regions in the phase diagrams where two phases with different symmetry of the 
superlattice (lamellar, hexagonal and bcc) coexist with each other in the finite temperature 
range. Since the block lengths in the system are assumed very large, this approximation 
becomes more accurate and we can ignore the fluctuation corrections |L2] . 
The paper is organized as follows: In section II, the model and the block length distribution 
are explained in details. Using the block length distribution one can calculate the correlation 
functions which are required for determination of the vertex functions of Landau free energy 
of multiblock copolymer. In section III we calculate the higher correlation functions up to 
sixth order and then the effects of the Laser distribution and Schultz distribution on them are 
compared. In section IV we give a general expression for the Landau free energy of a broad 
class of polydisperse multiblock copolymers. In Section V the influence of the parameters 
M and K on phase diagrams of AB multiblock copolymers are discussed. The paper ends 
with conclusion and two appendices. 



II Laser distribution and its polydispersity 

We consider an incompressible multiblock copolymer melts consisting of two types of blocks 
A and B in which length of both A and B blocks obeys Laser distribution. Generally the 
block length distribution is the statistical representation of the probability of attachment of 
one block to another block. The Laser distribution is used for the first time as a comprehen- 
sive block length distribution in multiblock copolymers in Reference , such that it can be 
reduced to other distributions in special cases. The normalized Laser distribution is defined 
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as 



0: 



, , K(M + n k-\ —K Kn(M + l) s ,2Kn?(M + 1)5 
P n = W — 2 exp iu 7 ^-i ^ H' 

where K > and is modified Bessel function. For K = 1, the Laser distribution reduces 
to the Laguerre distribution. As M — > oo, the Laser distribution reduces to Schultz-Zimm 
distribution [0, [0J. Finally for K = 1 and in the limit of M — > oo, the Laser distribution 
reduces to the Flory distribution. As it will be shown in Appendix I, the polydispersity of 
the Laser distribution is equal to kim+i) 2 • Clearly as M becomes infinite, its polydispersity 



reduces to j^. 



Ill Higher correlation functions 

In this section, we calculate the higher correlation functions, that is, the fourth, the fifth 
and the sixth order correlation functions. The calculation of the lower correlation function 
are given in details in Appendix II. The higher order correlation functions are defined as: 



,3, 3(1 -a) | f-2a' | a" 

y 4 y 3 2y 2 ' 



9aaaa = 4!(n A + n B Y( — + — — — + — ) 

9AAAAA = 5\(n A + n B 



4 .4(a-l) /-3a' a" a"\ 
(— + 1 ) 
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5 3(a-l) / + 4a' 3a" a'" a" 



9aaaaaa = 6!(n A + n B ) ( h — T - 5-r + 
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with 



— 1 r 

a® = ( Y / dnn 1 exp(—nx)P(n) 

Ha +riR J 



'n A + n B ' 

where a^ corresponds to a, a', a", a'", a"" for i=0, 1,2,3,4, respectively . As it is shown in the 
above equation, the value of a ^depends on the distribution function P(n) where here we use 
the most comprehensive distribution, that is, Laser distribution, where it can be reduced to 
Schultz distribution for paticular values of its parameters. In order to obtain the correlation 
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functions, we have to calculate the Laplace transforms of a(y). The Laplace transforms of 
a(y) (corresponding to Laser distribution) are: 

a® =RW.f3 k - 1 1 ~^Lt 1 \-— ) 

7 

k(M + l)*P , 1 ^ r k 
R = — — 7T — ( — ) 2 expl— ) 

_ fc(M + 1) + Mn n 
7 ~~ Mn n 

and 

/» = * (M+ , 1) , 

Mnl 

where the expression L,\ k is the Laguerre function. The influence of the parameters K 
and M on the shape of the second, the third and the sixth order correlation functions are 
given in the figures la to lc. The form of the third, fourth, fifth and sixth order correlations 
are same but their scales are very different( they are of order: 10 9 , 10 15 , 10 21 and 10 27 , 
respectively). If we apply the Schultz distribution, we obtain only the first point of the right 
hand of the curves. It should be reminded that as M — ► oo, the Laser distribution reduces 
to Schultz-Zimm distribution. 

IV Theory of Landau free energy 

From the point of view of Landau's phase transition theory ]rjj , a system under microphase 



separation undergoes a set of phase transitions which are described in terms of the order 
parameter: 

ip(x) = pa(x)~ < Pa >, 

where Pa{%) is the local density of A monomers and < pa > is the average value of the 
density over the volume V of the system. At low temperature, the order parameter can 
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become a spatially periodic function, which possesses the symmetries of a certain space 
group C i-e. : 

^(9) = S Q G LA? ex P (&0)8(q- Q), 

where ip{q) is the Fourier transform of concentration profile ij)(x), L is the reciprocal of the 
lattice L describing the translational symmetries of (, the amplitudes Aq and phases <3>q are 
adjustable parameters. The free energy of the system is : 

F(x,f) = -K B T\nZ(xJ), 

where T is temperature, K B is Boltzmann's constant, x is the Flory-Huggins parameter, f 
is volume fraction of the A block and Z is the partition function defined as the integral over 
all possible profiles ip(x), i.e. : 



di/>exp(-H(il>(x))/K B T), 

where the effective Hamiltonian, H(ip(x)) is the virtual free energy of the state with a given 
ip(x). Since we assume that Pa{x) differs from < p A > only slightly, the free energy for each 
symmetries of the superlattice i may be approximated by its Landau's expansion in powers 
of if)(x), i.e. : 

Fl(ip(x)) = K B T^2— [ dx x ■ •■dx n Y n (x 1 ■ ••x n )ip(xi) ■ ■■ip(x n ), 

with vertices T n defined as : 

r 2 (9i, 92) = V5 k {q u q 2 )[ 2 - 2 X , 

9AA{qi) * 9BB{qi) - 9ab{Qi) 

rs(gi , 02, qi) = -V5 k (qi, 93) 2^7(91 > 92, ^^(gi^feKfe), 

where g a/ 3 and are second and third correlation functions, respectively, and 

zM = 9a\{q) - 9aB(q)- 
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After inserting the trial function ifj(x) in Landau free energy, the lowest order term F 2 attains 
the form: 

This expression is independent of the structure. The third order and fourth order contribu- 
tion to free energy are different for the various structures: 

F£ ex 2 I 73 I A 3 

~ 3^ ' 
F| cc -4 I 73 I A 3 
~V~ ~ 3^6 ' 



while the fourth order contribution to the free energy has the form 

^ r 4(gi, 92, Qs, qt)A 4 . 

' 91,92,93,34 

In order to obtain the vertices of multiblock copolymer whose blocks lengths obey a certain 
distribution P(n), we have to calculate the correlation functions. The general expressions 
for the correlation functions are investigated in Appendix II. The phase diagram can be 
obtained by numerical minimization of free energy with respect to the parameters A and y. 
Several phases with different superlattice structure can coexist with each other. Considering 
a two-phase state with the volume fractions x and (1 — x) occupied by phases 1 and 2, 
respectively, the final expression for the free energy can be defined as : 

F = xF 1 (A 1 ,y 1 ) + (1 - x)F 2 (A 2 , y 2 ) - ^r(l - x)J 

and 

J = 4(— + l)A\ + 8(— ^— + \)A\A\ - 4(^ + \)A\. 

Here the phase diagram can be obtained by numerical minimization of free energy with 
respect to the parameters Ai, A 2 , yi and y 2 and x for each coexistence. 
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V The influence of the polydispersity on the phase di- 
agrams and phase coexistence 

In this section phase coexistence for various values of polydispersity with Laser distribution 
are presented in the framework of the mean-field theory. Here diagrams in the immediate 
vicinity of the critical point of the multiblock copolymers for K = 1 to 40 and M = 0.01 to 
100 are constructed. Figure 2 shows the mean-field phase diagrams for the most polydispers 
and the most monodisperse systems, respectively. In figure 2a, M is very large and K is 
very small and in figure 26, M is very small and K is very large. Figures 3 and 4 show 
the mean-field phase diagram for the multiblock copolymer melts in terms of variables N%, 
where x is the interaction parameter and f is the monomer fraction. In these figures system 
is at very high polydisperse and very low polydiperse states , respectively . It should be 
reminded that according to the equation of the polydispersity of Laser distribution, the 
polydispersity increases with increasing of M and decreasing of K. These figures show 
that the region of the stability of the phase coexistence shift considerably to the larger N% 
with increasing polydispersity. Figures 5 and 6 show the one-dimensional phase diagrams 
of multiblock copolymer with f = 0.25 for K = 1 and K = 40, respectively. One can see 
that for low values of K, the value of M has not considerable effect on the stability region of 
microphases while its effect is very large for higher values of K. Therefore, for large values 
of K, with the decreasing of M, the transition lines of the phases shifts severely upwards, 
while for small values of K these changes are very slow. Great attention should be paid to 
the scales in this figure. Even though these scales are unphysical, but it can be very useful 
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in showing the sensitivity of N% to M. The dependence of the N% on polydispersity, in the 
high polydisperse and low polydisperse systems, with f = 0.25, are given in figures 7 and 8, 
respectively. These figures show that for small values of K with increasing of M the slope 
of the curve is slow, while for large values of K, it is very sharp. 
CONCLUSION : 

In the present paper we have introduced the most comprehensive distribution function for 
the polydispersity of molecular weight distributions that is the Laser distribution. Laser 
distribution with two parameters K and M can be reduced to other distributions. The Flory 
distribution and Schultz-Zimm distribution are M — > oo limiting case of Laser distribution, 
therefore, so far the role of finite values of M has not been considered. Actually in this work 
we see that the effect of M (especially for the low M) on the curve of distribution versus block 
length and also its effect on the phase behaviour of multiblock copolymer is very interesting. 



Appendix I : The polydispersity of distributions 
Here in this appendix we derive the polydispersity of distribution functions. According to 
section III the polydispersity is defined as : 



U = 

where 



<n 2 > - <n> 2 



<n> 2 



< n* >= J dnP(n)n\ 



Using the above formulas we calculate the polydispersity of Laser distribution which is com- 
prehensive enough to include the other known distributions as a limiting cases. Hence we can 
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obtain the polydispersity of other distributions from this one. In calculating polydispersity 
of Laser distribution we have to use the following integral: 

roc , 

/ dxx p+v/2 exp(-ax)I v (2(3x*) = p\{3 v exp (p 2 / 'a)a { - p - v ~ 1) 'L p v (-(3 2 1 a) 
Jo 

. Using the above integral we obtain the following expression for the p-moment of laser 
distribution: 

<x p >=p\a- p L k p ~ 1 {-(3 2 /a), 

where 

L v p (x) = (l/p\)x~ v exp{x){d/dx) p {x v+p exp (-x)) 

is the Laguerre polynomial of order p. Now, using the first and second moments we obtain 
the following expression for polydispersity of Laser distribution: 

_ M(M + 2) 

K(M + iy 

Obviously for K = 1 we get the Laguerre distribution with polydispersity U = 1 — ( M +^2 
and as M ^ oo we obtain the Schultz-Zimm distribution with polydispersity U = where 
Flory distribution with polydispersity U — 1 can be obtained from it for K — 1. 



Appendix II: The Correlation Functions of Multiblock Copolymer [|T7j 
In a polydisperse multiblock copolymer melt several molecule types are presented, if s 
denotes a molecule type (i.e. a finite sequence of A-blocks and B-blocks) and if p s denotes 
the number density of molecules of type s, then the composition of the melt is fixed completely 
by the set p s . A molecule s is fixed if the lengths of all blocks in the molecule are fixed: 

8 = {NX,Nh,N%,Nl,...,N2,Nfj) 
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where N l a denotes the length of the i th h\ock of type a in the molecule. In this work we have 
assumed that the number of molecule (Q) is very large. Here P a {n) denotes the probability 
that a block of type a has a length between n and n+dn. For this distribution we have : 



/ dnP a {n) = 1 
Jo 



and 



roo 

/ dnP a (n) = n a , 
Jo 

where n a is the average length of a block of type a. The quantity p s can be defined as: 

Ps = cP s = cPa ( iV"^ ) P B ( N B ) Pa ( N\ ) P B ( N% ) . . . P A (N$)P B (N$). 

The normalization constant c is determined by the constraint that the monomer density is 
unity, i.e. : 

s 

where N s is the total molecule length. Combining the above equations gives the normaliza- 
tion constant c: 



Q(n A + n B )' 

Now we calculate the second order function. The average second order function over the 
composition p s is given by: 

9ap = ^p s ^a i (a)a i ( / 9)exp(-g 2 \ j - i \) 

s i,j 

where <Ti(a) equals 1 only if block i is of type a and otherwise is zero. Where i and j are 
many blocks apart, the weight factor exp(— q 2 \ j — i |) will give a negligible contribution to 
the correlation function, and one gets a contribution only if the number of blocks between i 
and j is of order unity. Here for simplicity, the calculation of the correlation function where 
i and j are in the same block are given. For second correlation function, the AA-component 
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of the second order correlation function is equal to : 



9aa = 



Q(n A + n B ) 



/rn rn 
dnP(n) di dj exp(— x(j — i)) — 
Jo Jo 

, N n exp(— nx) 1 . , . — 2(1 — a(y)) 2f . 

dnP(n)- + y \ J - — = (n A + n B )( K KU " + 
ua + n,B J x x z x z y z y 

where: 

f— ,x = q 2 ,y = (tia + ub)x, a(y) = dnexp(—nx)P(n) 

WA + Ub J 

For the BB-component of the second order correlation function, j3(y) is used to as the 
Laplace transform of the length distribution of the B-blocks. Since the states corresponding 
to i < j give the same contribution as j < i ones, one should take into account a factor 2 
due to this symmetry. Finally the factor Q can be interpreted as the number of blocks in 
which monomer i can be presented. The contribution third order correlation function for the 
situation where all three monomer are in same block A is : 

| roo rn rn rn 

9AAA = -pr, ; r3!Q / dnP(n) di dj dkexp(-x(k - i)) 

Q(nA + nB) Jo Jo Ji Jj 



at , \2/ 2(1 -a) /-a, 
Q(n A + n B ) (- 



y 3 y 2 



a is the derivative of a with respect to y: 

roo 

a{y) = d/dy \ dnexp(—nx)P(n) 
Jo 

dnnexp(—nx)P(n) 



io 

oo 



(n A + n B ) Jo 
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Figures captions 

FIG. 1. Dependence of the second, third and sixth order correlation functions on the value 
of K for various values of M, respectively (la-lc). 

FIG. 2. Mean-field phase diagrams in the vicinity of the critical point for multiblock copoly- 
mers for the most polydispers and the most monodisperse system, respectively. The lower 
line separates the disordered phase from the bcc phase, the middle line separates the bcc 
phase from the hexagonal phase and the upper line separates the hexagonal phase from the 
lamellar phase. (2a) K = 1 and M = 100. (26) K = 40 and M = 0.01. 
FIG. 3. Mean-field phase diagram in the vicinity of the critical point for multiblock copoly- 
mers with K = 1 and M = 100. White areas denote single phase regions lam, hex and bcc 
from up to down, respectively. Black areas denote the regions of coexistence of 2-phases. 
FIG. 4. Mean-field phase diagram in the vicinity of the critical point for multiblock copoly- 
mers with K = 40 and M = 0.01. White areas denote single phase regions lam, hex and 
bcc from up to down, respectively. Black areas denote the regions of coexistence of 2-phases. 
FIG. 5. The one-dimensional phase diagrams of the multiblock copolymer melts f = 0.25 
for the Laser distribution with K = 1 for various values of M. 

FIG. 6. The one-dimensional phase diagrams of the multiblock copolymer melts f = 0.25 
for the Laser distribution with K = 40 for various values of M. 

FIG. 7. Dependence of the parameter N% on the polydispersity from M = 0.01 to M — ► oo 
for K = 1. 

FIG. 8. Dependence of the parameter Nx on the polydispersity from M = 0.01 to M — > oo 
for K = 40. 



